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Radiative multipole moments of scalar, electromagnetic, and linearized gravitational fields in 
Schwarzschild spacetime are computed to third order in v in a weak-field, slow-motion approxima- 
tion, where v is a characteristic velocity associated with the motion of the source. These moments are 
defined for all three types of radiation by relations of the form ^(t, x) = r~ ^2 lm Mim(u) Yi m (9, <f>), 
where $ is the radiation field at infinity and Mi m are the radiative moments, functions of retarded 
time u = t — r — 2Mln(r/2M — 1); M is the mass parameter of the Schwarzschild spacetime and 
(t,x) = (t,r,9,<f>) are the usual Schwarzschild coordinates. For all three types of radiation the 
moments share the same mathematical structure: 

To zeroth order in v, the radiative moments are given by relations of the form Mi m {u) oc 
(d/du) 1 J p(u,x)r l Yi m (d,<t>) dx, where p is the source of the radiation. A radiative moment of 
order / is therefore given by the corresponding source moment differentiated I times with respect to 
retarded time. 

To second order in v, additional terms appear inside the spatial integrals, and the radiative 
moments become M.i m (u) oc (d/du) 1 J [1 + 0(r 2 d 2 ) + O(Mfr)] pr l Yi m dx. The term involving 
r 2 d 2 can be interpreted as a special-relativistic correction to the wave-generation problem. The 
term involving M/r comes from general relativity. These correction terms of order v 2 are near- 
zone corrections which depend on the detailed behavior of the source. Furthermore, the radiative 
multipole moments are still local functions of u, as they depend on the state of the source at retarded 
time u only. 

To third order in v, the radiative moments become Mi m {u) — > Mi m (u) + 2M J_ [ln(w — u) + 
const] wMj m (ii') du , where dots indicate differentiation with respect to u . This expression shows 
that the 0(v 3 ) correction terms occur outside the spatial integrals, so that they do not depend on 
the detailed behavior of the source. Furthermore, the radiative multipole moments now display a 
nonlocality in time, as they depend on the state of the source at all times prior to the retarded time 
u, with the factor \n(u — u) assigning most of the weight to the source's recent past. (The term 
involving the constant is actually local.) The correction terms of order u 3 are wave-propagation 
corrections which are heuristically understood as arising from the scattering of the radiation by the 
spacetime curvature surrounding the source. 

The radiative multipole moments are computed explicitly for all three types of radiation by 
taking advantage of the symmetries of the Schwarzschild metric to separate the variables in the 
wave equations. Our calculations show that the truly nonlocal wave-propagation correction — the 
term involving \n(u — u) — takes a universal form which is independent of multipole order and field 
type. We also show that in general relativity, temporal and spatial curvatures contribute equally 
to the wave-propagation corrections. Finally, we produce an alternative derivation of the radiative 
moments of a scalar field based on the retarded Green's function of DeWitt and Brehme. This 
calculation shows that the tail part of the Green's function is entirely responsible for the wave- 
propagation corrections in the radiative moments. 

PACS numbers: 04.25.Nx, 04.30.-w, 04.30.Db, 04.40.-b 



I. INTRODUCTION AND SUMMARY 



A. Tails in waves 



It has long been known that in general, the propagation of massless fields in curved spacetime does not proceed 
along characteristics only, but is accompanied by wave tails. Much attention has been devoted to this topic since the 
ground-breaking work by Hadamard [Q. Here are some of the highlights: 

In 1952, Choquet-Bruhat Q studied the initial value problem of general relativity and showed that the gravitational 
field at some event P depends not only on the data put on the intersection of P's past light cone with the initial 
surface, but also on the data put inside this region. This result indicates that in general relativity, field propagation 
proceeds at all speeds less than, or equal to, the speed of light. 
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In 1960, DeWitt and Brehme || constructed Green's functions for the scalar and electromagnetic wave equations 
in arbitrary curved spacetimes, and showed that these split naturally into a direct part, with support on, and only on, 
the light cone, and a tail part, with support within the light cone. A similar Green's function was also constructed 
for the Einstein field equations 0. 

In 1968, Kundt and Newman [p| established that for hyperbolic partial differential equations in two dimensions, 
the presence of wave tails is the rule rather than the exception. This conclusion was extended by McLenaghan and 
co-workers to the case of conformally invariant wave equations in four dimensions. 

Wave tails are known to have important physical consequences. For example, DeWitt and Brehme || have shown 
that the tail part of the electromagnetic field is of paramount importance in deriving the equations of motion for 
charged particles in curved spacetime. Similarly, Mino, Sasaki, and Tanaka [lCf| , as well as Quinn and Wald [ylj , 
have recently shown that tails are entirely responsible for the gravitational radiation reaction force. And as a final 
example, Price fll2| has shown that gravitational-wave tails play an integral part in the physical process by which a 
recently formed black hole relaxes to a stationary state, as is demanded by the no-hair theorems. 

The presence of tails in the gravitational waves produced by an isolated source was first demonstrated in 1965 
by Bonnor and Rotenberg jyj). In 1968, this work was extended by Couch, Torrence, Janis, and Newman who 
showed that an initially outgoing wave will be partly backscattered by the spacetime curvature surrounding the source, 
thereby creating an incoming wave. A further extension of this work has appeared very recently |l5f| . 

In 1992, Blanchet and Damour |l(| considered for the first time the effect of tails on the behavior of the gravitational 
waves at infinity, thereby concentrating on effects that could potentially be observed directly. They found that the 
gravitational waves at time t depend not only on the state of the source at the corresponding retarded time u 
[essentially u = t — r/c, where r is the distance to the source, but see Eq. (1.2) for a more precise definition], but also 
on the state of the source at all times prior to the retarded time. (Once again, this indicates that wave propagation 
proceeds at all speeds less than, or equal to, the speed of light.) Subsequently |l7|-[20|], it was shown that tails play 
an important role in the generation of gravitational waves by the orbital inspiral of a compact binary system. These 
waves are among the most promising for detection by future kilometer-scale interferometers such as the American 
LIGO (Laser Interferometer Gravitational- wave Observatory J2l|]) and the French- Italian VIRGO 

Given the physical relevance of tails in the propagation of radiation in curved spacetime, it appeared to us worthwhile 
to seek a deeper understanding of this phenomenon by asking how it depends on the type of radiation being considered, 
and by digging further into the nature of its physical origin. This paper reports on the results of such an investigation, 
in which we study the influence of tails on those properties of massless fields that are directly measurable to an observer 
at infinity: the radiative multipole moments. We consider the cases of scalar, electromagnetic, and gravitational 
radiation generated by an isolated source, and propagating to infinity in a spacetime curved by a nonrotating central 
mass M . All of our results are derived on the basis of a weak- field, slow-motion approximation. Throughout this 
paper we use units such that G = c = 1, and we employ the definitions and conventions of Misner, Thorne, and 
Wheeler M. 



B. Scalar radiation 



We begin the summary of our results with the simplest case, that of a scalar field <!>(a:) obeying the wave equation 
□ $ = — 47rp, where □ = g a ^V a V p, g a p is the Schwarzschild metric, and p(x) a given source. The symbol x 
collectively designates all Schwarzschild coordinates {t, r, 8, <fi}. As is shown in Sec. Ill, the radiative part of the scalar 
field, which dominates at infinity, can be written as 

oo / 

$rad(*,z) = ~^ Y, Zlm(u)Y lm {e,<t>), (1.1) 
2=0 m=-l 

where Zi m (u) are radiative multipole moments, depending on retarded time 

u = t-r-2Mln(r/2M-l), (1.2) 

and Yi m (9, (f) are the usual spherical harmonics. The symbol x collectively designates all spatial coordinates {r, 6, </>}. 

In a leading-order calculation in a weak-field, slow-motion approximation, the radiative multipole moments are 
found to be given by 

(azTl)n(^) / P(«>ViWM)d*, (1.3) 
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where dx = r 2 dr d cos 9 d<fi and the int egration is over the region of space occupied by the source; we assume that 
this region is bounded. Equation (E3) shows that the radiative moments are obtained from the source moments 



/ pr l Y[ m dx by taking a number of time derivatives equal to the multipole order. 

In a more accurate calculation, incorporating corrections of order v 2 with respect to the leading-order results (with 
d« 1 a characteristic velocity associated with the motion of the source), we find 



4?r 



(2Z + 1)!! \du 



1 



(rd u ) 



,M 



2(2^ + 3) 



l- 



p[u, x) r l Y lri 



I dx. 



(1.4) 



It is easy to see that the correction terms are indeed of order v 2 . The term involving [rdu) is of order (r c /t c ) 2 , where 
r c is a characteristic radius within the source and t c a characteristic time scale associated with its motion; the ratio 
r c /t c defines the characteristic velocity v. This term can be understood as arising from special-relativistic corrections 
to the wave-generation problem. On the other hand, the term proportional to M/r comes from general relativity, and 
is also of order v 2 by virtue of the virial theorem for bound motion in a gravitational field. It should be noted that 
in Eq. (1.4), the correction terms occur inside the spatial integrals, so that they depend on the detailed behavior of 
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the source. Furthermore, these corrections are purely local in time, as (u) depends on the state of the source at 
the time u only. The terms of order v 2 are therefore near-zone corrections that have nothing to do with the tail effect 
discussed previously. 

A calculation carried out to order v 3 in a weak-field, slow-motion approximation does reveal the influence of the 
tails. Indeed, the radiative multipole moments are now given by 



Z^[u)=Z^[u)+2M 



(2), 



In 



AM 



scalar 



+ 7 



which clearly displays a nonlocality in time. Here, a dot indicates differentiation with respect to v! , and 



^scalar = ^ + y + 



2' 



(1.5) 



(1.6) 



where tp[l + 1) = —7 + J2k=i ^ s ^ ne digamma function (7 ~ 0.57721 is Euler's number). We notice that the 
correction terms of order v 3 occur outside the spatial integrals, so that they do not depend on the detailed behavior 
of the source. These are wave-propagation corrections, which are readily associa ted with the occurrence of wave tails 
in curved spacetime. We must point out that of the three 0(v 3 ) terms in Eq. 

Jscalar i 
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two are actually instantaneous, 

as they are equal to 2M(/3 ; scalar + 7)2^ (it). The remaining term is truly nonlocal, and the factor ln(u — u') assigns 
most of the weigh t to the source's recent past. 

Equations (pj)-(p^) were derived by integrating the wave equation for a scalar field in Schwarzschild spacetime, for 
which the variables can conveniently be separated. (The actual work of integrating the radial equation is carried out 
in Sec. II.) But because our calculations use only the weak-field behavior of the Schwarzschild solution, our results are 
insensitive to the detailed form of the metric. Although we rely heavily on the symmetries of the Schwarzschild metric 
to separate the variables in the wave equation, our results rely only on the fact that the field is spherically symmetric 
at large distances. Our results would therefore hold also in more general spacetimes, with spherical symmetry holding 
only approximately at large distances. Staticity, however, is a crucial assumption, and our results would not be valid 
if the spacetime were rotating. Alt hough the wave-propagation corrections to the radiative multipole moments would 
take the same form as in Eq. (1.5), the spacetime's rotation would create additional terms of order v 3 . These would 



occur inside the spatial integrals, and would describe near-zone corrections of the spin-orbit type 



C. Electromagnetic radiation 

In Sec. IV we turn to the case of electromagnetic radiation produced by a given current density J" (a;) in 
Schwarzschild spacetime. (The remarks of the preceding paragraph, regarding the generality of our results, apply 
equally well here.) The radiative part of the vector potential is given by 

A ^ = - E E Y«' lm & *) + Slm[u) Y*> lm [e, 0)] , (1.7) 

1=1 rn=—l 

where Xi m (u) and Si m [u) are charge and current multipole moments, respectively, while Y^' lm (9, if) and Y^' lm (9, tf>) 
are the vectorial spherical harmonics described in Appendix A. In a calculation accurate to order v 3 in a weak-field, 
slow-motion approximation, we find 
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In 



4M 



+ A em + 7 



(1.8) 



and 



S%(u)=S%(u)+2M 



In 



4A/ 



+ /?f m + 7 



(1.9) 



where 



/3f m = Ml + 1) + 



1 



2 21(1 + 1) 



(1.10) 



These relations are analogous to Eq. fll.q) and have the same physical interpretation. The second-order multipolc 
moments are given by 



*2(«) 



4tt 



Z + 1 / d 



(2/ + 1)!! V Z V du 



(r(u,f)r'ij m (0,0)di? (1.11) 



and 



S lm( u ) 



4?r 



(2Z + 1)!! \du 



(rd u ) 



P - 1 M 



2(2Z + 3) Z r 



J a {u,xyY*> im {6,(j))dx 



(1.12) 



where a = J f — rd u J r /(I + 1). These relations are analogous to Eq. (1.4) and have the same physical interpretation. 
Finally, the zeroth-order moments, 2^(u) and S^(u), are obtained from Eqs. ( 1.11 ) and ( 1.12 ) by discarding all 
0(v 2 ) terms; in this limit, a = J*. 



D. Gravitational radiation 



The case of gravitational radiation is conceptually very different from the previous cases, because of the fact that the 
spacetime metric is now dynamical. However, if we assume that T a "(x), the stress-energy tensor responsible for the 
radiation, is small, then the Einstein field equations may be linearized in the small deviation of the metric with respect 
to the Schwarzschild form. This results in a wave equation for the metric perturbation [ p7| , and mathematically, the 
gravitational-radiation problem ends up resembling closely the scalar and electromagnetic analogues. This is the 
problem considered in Sec. V. 

The traceless-transverse gravitational- wave field is given by |^8| 



1=2 m=-l 



,(u) T™ M (0, 0) + S lm (u) T™< lm (0, 



(1.13) 



where Xi m (u) and Si m (u) are mass and current multipole moments, respectively, while T a *' m (6, (f) and T^' lm (9, eft) 
are the tensorial spherical harmonics described in Appendix A. The weak-field, slow-motion approximations to the 
multipole moments are 



. 11 U \ /^grav.mass 



2jm («')**' 



(1.14) 



and 



S%( u )=S%(u) + 2M 



.11 11 \ ^grav, current , 

ln, ^iH +/3f +7 



(1.15) 



where 



.grav.cuxrent = ^ + ^ + 



2 1(1 + 1) 



(1.16) 



and 
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rav,mass ^grav, current 



(i-l)Z(J + !)(! + 2)' 



(1.17) 



Equations ( 1.14 )-( 1.16| ) were previously derived (in a different representation involving symmetric tracefree tensors) 
by Blanchet ]29fl (see also Refs. and [fL9[ ) on the basis of post-Newtonian theory. Blanchet, however, derived these 
relations without an explicit knowledge of the second-order moments, which are given by 



*S(«) 



16tt 



'(/ + !)(/ 



(2Z + 1)!! V 2(1 - 1)1 



2(1- 



1+1 (rd ) 2 -(l + 2)^ 

-l)(2Z + 3) 1 u) [l + Z) r 



a(u, x) r l Yi r< 



\dx (1.18) 



and 



5g(«) 



16tt 



(21 + 1)!! 



2(Z + 2) 
Z - 1 



(Z- l)(Z + 2) M 
Z V 



1 + 4 



2(Z + 2)(2Z + 3) 



J a (u,f)r'y c 



l vB.lr, 



(1.19) 



where a = T* 



r 2 T° 



9 • 2 

r sin 



- Ard u T tr /(l + 1), J a = T ta - rd u T ra /(l + 2). To our knowledge, 
such explicit expressions for the radiative multipole moments have not appeared before in the literature. Finally, the 
zeroth-order expressions are recovered by discarding all 0(v 2 ) terms from Eqs. (LIS) and ( l.lSj ); in this limit, a = T u 
and J a = T ta . 

The physical interpretation of these results is the same as in the previous cases, and the results share the same 
degree of generality as the previous ones (see the concluding paragraph of Sec. I B). 



E. Universality of the tail correction 



A survey of the preceding subsections reveals that the multipole moments of scalar, electromagnetic, and grav- 
itational radiation fields all share the same mathematical structure, with terms of order v 2 near-zone corrections 
depending on the detailed behavior of the source, and with terms of order v 3 wave-propagation corrections indepen- 
dent of the detailed behavior of the source. And while the near-zone corrections are local in time, the wave-propagation 
corrections are not. 

We also observe that the tail corrections depend on the multipole order I and on the field's type only through the 
terms involving the various /Vs. These terms are actually instantaneous, because after integration over du' ', they are 
found to be proportional to the first derivative of the zeroth-order moments evaluated at u. The truly nonlocal tail 
corrections, which involve the weighting factor ln(u — u'), are independent of multipole order and field type. This 
remarkable result, that the tail correction has a universal form, is one of the main new contributions of this paper. 



F. Physical origin of the tail term 



The nonlocality (in time) of the radiative multipole moments is heuristically understood as arising from the scat- 
tering of the radiation field by the spacetime curvature surrounding the mass M , and a survey of our previous results 
reveals that the tail terms are indeed proportional to M. Now, the mass parameter enters twice in the metric of an 
asymptotically flat spacetime: Assuming that the weak-field metric is expressed in Schwarzschild-like coordinates, we 
have g u ~ — 1 + 2M/r and g rr ~ 1 + 2M/r at large distances. Because of this degeneracy, it is impossible to tell 
whether it is "gtt$ mass" which is "mostly responsible" for the tail effect, or whether it is "g rr 's mass", or whether 
both are "equally responsible". In other words, we cannot tell how the temporal and spatial curvatures separately 
contribute to the tail effect. 

We examine this question in Sec. VI, for the specific case of scalar radiation. To lift the degeneracy, we artificially 
introduce an additional mass parameter, 7M, in the description of our spacetime. This is defined so that the metric 
functions are now given by gu ~ —1 + 2M/r and g rr ~ 1 + 2-yM/r in the weak-field limit. General relativity is 
recovered by setting 7 = 1. 



Integrating the scalar wave equation for the modified spacetime yields Eq. (LI) with 



( 3 )f»A - srCO 



Z^(u) + (l + l)M 



In 



scala 



+ 7 



(1.20) 
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where f3f 



scalar 



^+1)+ 7/(1+7), 



4™» 



4n 



(2/ + 1)1! 



(rd u ) 2 (2/ + 1) 7 -1 M' 



2(2? + 3) 



p(u,x)r l Yi m {6,4>)dx, 



(1.21) 



and Z^(u) is given by Eq. (1.3). These expressions reduce to Eqs. ( |L5| ) and (1.4), respectively, when 7 = 1. We 
see that in the modified spacetime, the tail terms are proportional to (1 + 7). This allows us to conclude that in 
general relativity, the temporal and spatial curvatures contribute equally to the tail effect. This intriguing observation 
is another main contribution of this paper. 



G. Spacetime approach 



The final section of the paper, Sec. VII, contains an alternative derivation of Eqs. (1.1 ) — (pi . 5|) based on the spacetime 
approach of DeWitt and Brehme [0. For simplicity, we again restrict attention to the case of scalar radiation. 

The mathematical methods employed in Sees. Ill, IV, and V to derive expressions for the radiative multipole 
moments of integer-spin fields are based upon a separation of variables approach made possible by the symmetries of 
the Schwarzschild solution. These methods, though convenient for practical computations, do not reflect closely the 
physical picture of wave propagation in curved spacetime. In particular, the distinction between direct terms (which 
are local in time) and tail terms (which are nonlocal) emerges only at the very end of the calculation. 

The spacetime approach of Sec. VII is based instead on G(x,x'), the DeWitt-Brehme retarded Green's function 
for the scalar wave equation Q. As was mentioned in Sec. I A, G(x, x') is expressed as a sum of two parts. The first 
part has support on, and only on, the past light cone of the field point x, and gives rise to the direct terms in the 
waves. The second part has support inside the past light cone of x, and gives rise to the tail terms. In the spacetime 
approach of DeWitt and Brehme, the mathematics reflect the physical picture quite closely. 

The radiative multipole moments calculated in Sec. VII agree precisely with those calculated in Sec. III. Therefore, 
the calculation based on the spacetime approach tells us nothing new in terms of the final answer. Nevertheless, this 
alternative derivation is very instructive, because of the fact that the mathematical origin of the tail correction is 
clear from the outset. We regard this as another important contribution of this paper. 



H. Organization of this paper 



The remaining sections of the paper contain the detailed derivations of the results summarized above. After 
laying some preliminary ground work in Sec. II, we integrate the wave equations for scalar, electromagnetic, and 
gravitational radiation in Schwarzschild spacetime in Sees. Ill, IV, and V, respectively. All calculations are carried 
out in a weak-field, slow-motion approximation. In Sec. VI we integrate the scalar wave equation for the artificially 
modified spacetime. And finally, in Sec. VII we integrate the scalar wave equation using the spacetime approach of 
DeWitt and Brehme S. Various technical details are relegated to five Appendices. 



II. GENERALIZED REGGE- WHEELER EQUATION 



The generalized Regge- Wheeler equation 



dr 



1(1 + 1) 2(s 2 -l)M 



>Xi(w;r) = 0, 



(2.1) 



has long been known to govern the evolution of integer-spin fields in Schwarzschild spacetime. Here, r is the usual 
Schwarzschild coordinate, / = 1 — 2M/r (with M denoting the mass of the spacetime), and d/dr* = fd/dr. Also, ui 
denotes the frequency of the field, I its multipole order, and s = {0, 1, 2} its spin. The precise relation between the 
mode functions Xi(u),r) and the corresponding scalar, electromagnetic, and gravitational fields will be described in 
Sees. Ill, IV, and V, respectively. In this section we consider the purely mathematical problem of integrating Eq. (2.1) 
in the low- frequency limit. 

We first examine the question of boundary conditions. It is easy to check that Xi(uj;t) must behave as e ±lujr , 
where 



G 



r* = r + 2Mhx(r/2M - 1), 



(2.2) 



in the asymptotic limits r — > 2Af, r — > oo. It will become clear in the following sections that the desired solution 
is the one which describes purely incoming waves at the black-hole event horizon. We therefore select the function 
Xf{uj] r), such that 



X"(lo; r -> 2M) ~ (const) e 



(2.3) 



The constant appearing in front of e~ lulr determines the overall normalization of the Regge- Wheeler function. Because 
our final results will be independent of this normalization, we shall leave this constant arbitrary. At infinity, X, (oj; r) 
describes a superposition of incoming and outgoing waves. Consequently, 



X" ( W ; r -> co) ~ AT(uj)e^ r + A^^Y' 



(2.4) 



The amplitudes A\ u (u>) and A? ut (uj) are determined by solving the differential equation. 

We wish to integrate Eq. (^J) in the low- frequency limit, for M\w\ <C 1. Without loss of generality, we henceforth 
take u> to be positive; the negative-frequency case can easily be recovered from the relation Xp (— u>) = X^(u>), where 
a bar denotes complex conjugation. To facilitate the calculations, we define the small (positive) quantity 

e = 2Mw, 

and introduce a new dependent variable, 

z = u>r. 

After substitution and expansion in powers of e, the generalized Regge- Wheeler equation becomes 



(2.5) 



(2.6) 




1(1 + 1) 



1 d 

z dz 



1(1 + 1) + 1 - s 2 



+ O(e 2 )}X l (z) = 0. 



(2.7) 



It should be noted that when expanding in powers of e, it was implicitly assumed that z 3> s . O ur low-frequency 



solution will therefore be restricted to the domain r ^> 2M. Our task is now to integrate Eq. fl2.7|) to first order in 



X^ + eX^ 1 ' + 0(e 2 ), substituting into Eq. (2/7), and solving order by 



e. We proceed by iteration, by writing Xi 
order. 

Because we are restricted to the domain z 3> e, which excludes the event horizon at z = e, Eq. (2.3) cannot be 
imposed directly, and the issue of boundary conditions must be re-examined. This question was addressed by Poisson 
and Sasaki |3l]] , who integrated the s = 2 Regge- Wheeler equation in the domain z 1 (which includes the horizon) , 
impose d th e correct boundary condition at z = e, and then matched the resulting function to the general solution 
of Eq. (2.7) in the common domain e < z < 1. Such an analysis will not be repeated here. It suffices to state the 
conclusion: To be compatible with the incoming- wave boundary condition at the horizon, the solution to Eq. (2.7) 
must be regular in the (unphysical, and unrealized) limit z 



0. 



With this in mind, the desired zeroth-order solution to Eq. (2.7) is 



X 



ff(0) 



(z) = zji(z), 



(2.8) 



where ji(z) are the spherical Bessel functions of the first kind. It should be noted that Eq. ( |2.g| ) provides a particular 
choice for the overall normalization of Xf 1 (z). The first-order solution is then determined by solving 



dz 2 



1(1 + 1) 



-Wi(z), 



where 



Wt(z) 



1 d 

z dz 



1(1 + 1) 



2 I 



zji(z)- 



The general solution to Eq. (2.£) is 

X\ 1 \z) = z jl (z) 



z'ni(z')Wi(z')dz' 



+ zni(z) 



z'3i(z')Wi(z')dz' 



(2.9) 



(2.10) 



(2.11) 
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where ni(z) are the spherical Bessel functions of the second kind, and a and b are constants which must be chosen so 
that X^\z) satisfies the r egula rity condition at z = 0. 

The integrations of Eq. ( 2.11 ) can be carried out explicitly. First, we use the recurrence relations among spherical 
Bessel functions (Ref. |32], p. 439) to write Wi(z) in the form 



wn 9 . , s (i-s)(i + s) . (i- s + m + s + i) . 

zWi(z) = 2zji{z) —— Ji-x(z) — — ■ 3l+i{z) 



21 + 1 



21 + 1 



Second, we evaluate the integrals using the results found in the Appendix of Ref. 
ulations, we arrive at 



(2.12) 

After straightforward manip- 



X\ 1] (z) =\a-A l (z)} zji (z) + \b + 7 + B l (z)] zn, (z) - ^ffj + g) zji-^z) 



(I- s + l)(l + s + l) 
2(i + l)(2i + l) 



zji+i(z), 



21(21 + 1) 



(2.13) 



where 



l-i , 

Mz) = Si(2z) + z 2 n (z)j {z) + ^ - 



~{\P P+l 



z 2 n p {z)j p (z) 



(2.14) 



and 



B t (z) = Ci(2z) - 7 - ln(2z) + z 2 j 2 (z) + ^ - 



^-f VP P+l 

p— l v 



Viz)- 



(2.15) 



Here, Si and Ci are the sine and cosine integral functions, respectively, and 7 ~ 0.57721 is Euler's number. 
In Appendix B, the functions Ai(z) and Bi(z) are evaluated in the limit z — > 0. We find 



Mz) = ~ + 0(z% 



Bi(z) = - 



z 



21+2 



+ 0{z 2l+i ) 



l(l + l)(2l-l)U(2l + iy.\ 



(2.16) 



It follows that lW(«) goes to zero in the limit z — > provided that = —7. Otherwise, the function diverges. This 
choice for 6 therefore selects X ; 1 (z), the desired solution. The constant a remains arbitrary, because it affects only 
the overall normalization of the solution; a can be set to zero without loss of generality. 

Integration of the generalized Regge- Wheeler function, to first order in e, is now completed. We have pointed out 
that our answer incorporates a specific choice of overall normalization which is provided by Eq. ( |2.8| ) . We now wish to 
form a normalization-independent quantity, X^(z)/A\ n , which will be required in the following sections of this paper. 
We must therefore calculate A\ n . This involves the evaluation of X^(z), a s giv en by Eqs. ( |2.8| ) and ( |2.13j )-(2.15), in 
the limit z — > 00, and a comparison with the low-frequency limit of Eq. (2.4). Such a calculation is carried out in 
Appendix C. The final result is quoted here: 



A) n 



= 2(-t) l+1 e^ ln2 "-^(l + fe)*| [l - eA l {zj\j l [z) + eBi(z)m(z) 

(l-s)(l + s) . . . (I - s + l)(l + s + 1) . , , 2 , 



21(21 + 1) 



2(Z + 1)(2Z + 1) 



(2.17) 



where 



A = W + 1) + 



1 



21(1 + 1)- 



(2.18) 



with t/j(l + 1) = —7 + k 1 denoting the digamma function. 



Finally, we evaluate Eq. ( |2.17 ) in the limit z <C 1, to a degree of accuracy sufficient for our purposes in the following 
sections of this paper. For this, we use Eq. (2.16) and the series expansions for the spherical Bessel functions. We 
arrive at 
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^(w) (2/ + 1) 



e 2iM w (ln4M|c|-/3 i ) ^ + ^M^Q (_ lwr ) 



1 + 1 



X U 



(wr) 2 + 



2(2/ + 3) J ) 

This equation holds both for positive and negative frequencies. 







(wr) 4 ,Mw 2 r,(M/r) s 



(2.19) 



III. SCALAR RADIATION 



A. Wave equation 

We begin our study of radiative multipole moments in curved spacetime with the simplest case, that of a real scalar 
field <5(cc) obeying the wave equation 



= -4wp(x). 



(3.1) 



Here, □ = g a ^W a W p is the curved spacetime wave operator, p(x) is an unspecified source function, and x collectively 
designates all spacetime coordinates. The spacetime is assumed to be Schwarzschild (with mass M), and the usual 
coordinates {i, r, 6*, </>} are adopted. 

Because the spacetime is static and spherically symmetric, the scalar field can be decomposed according to 



$(x) = ^J duj^2Ri m (uj;r) Y lm (d,< 



(3.2) 



where the sums over I and m are restricted by I > 0, |m| < I. Substituting this into Eq. (3.1), we obtain the following 
ordinary differential equation for the radial function Ri m (uj; r): 



dr* 



1(1 + 1) + 2M 



>Ri m (u>; r) = fTi m (u\ r), 



where d/dr* = fd/dr and / = 1 — 2M/r. The source term is given by 

Ti m (uj;r) = -47rr / p(u> , x) Yi m (0 , </)) dSl, 



where dft — d cos 6 dip, a bar denotes complex conjugation, and 

p(uj,x) = j- J p(t,x)e^dt 
is the Fourier transform of p(x). The symbol x collectively designates all spatial coordinates. 



(3.3) 



(3.4) 



(3.5) 



B. Solution 



Equation ( |3.3| ) has a Sturm-Liouville form, and it can therefore be solved in terms of a Green's function constructed 
from two linearly independent solutions to the homogeneous problem. Which solutions are selected depends on 
the boundary conditions we wish to impose on Ri m (u>;r). The appropriate choice here is dictated by the physical 
requirement that the scalar field must represent waves which are pur ely ingoing at the black hole horizon (r = 2M), 
and purely outgoing at r — oo. This amounts to integrating Eq. (3^) with a no- incoming-radiation initial condition. 

We therefore seek functions i?f f (w;r) and Rf°(uj;r), solutions to 



dr* 



1(1 + 1) 2M 



>Ri(u;r)=0, 



(3.6) 



and such that 







R?{us;r^ oo) ~ Q^Me"^* + Q{ ut {u)e™ r * , 
Rf°(u;;r -> oo) ~ e^ r *, 



(3.7) 



where r* = r + 2M ln(r /2M — 1). Equation (3.7) indicates that i? ; ff (w;r) describes waves which are purely ingoing 
at the black-hole horizon, while Rf°(uj;r) describes waves which are purely outgoing at infinity. The behavior of 
Rf°(ui; r) near r = 2M will not be needed. 



In terms of these functions, the solution to Eq. (3.3) takes the form 

1 



Rim(u;r) 



Rf°{u;r) I Ti m (u>;r')Rf(cj;r')dr' 

2M 



+ Rf(uj;r) I T lm (u-S)Rr > (w-S)dr' 



(3.8) 



where the factor 2iujQ\ n (u>) is the conserved Wronskian of the functions R?(u; r) and R^°(uj; r). 

We are interested in the radiative part of the field, which dominates at large distances from the source. More 
precisely, we define the radiative field by r<j> ra d(rc ) = hm r _ >0o r$>(x), which expresses the fact that at large distances, 
$(cc) = $ ra( j(a;) + 0(l/r 2 ). Evaluating Eq. ( [3.8[ ) in the limit r — > oo [assuming that Tj m (w;r) has compact support], 
we obtain 



Rim{u m , r —y oo) ~ Zi m (u>) e l 



where 



2iwQ}?(u) J 2M 
Finally, substituting this into Eq. ( |3.2|) yields 

$rad(t, x) = ~Yl Z lm {u)Yi 



Ti m (u;r)Rf(uj;r) dr. 



Irn \ 



! in 



where u = t — r* is retarded time, and 



Z lm {u) = / Z lm (iu)e' MU dco. 



(3.9) 



(3.10) 



(3.H) 



(3.12) 



The quantities Z\ m (u), or their Fourier transform Zi m (uj), will be referred to as the radiative multipole moments of 
the scalar field $>(x). 



C. Slow-motion approximation 



In Eq. (pTToj) , the radiative multipole moments are written in exact form in terms of the source T/ m (aj; r) and the 
function i?7 r (w;r)/Qj n (w). While the source function will be left unspecified, we now wish to find an expression for 



Rf 1 (uj; r)/Q\ n (uj). To do this we must resort to approximations, because Eq. (3.6) cannot be integrated in closed form. 

We will derive approximate expressions for the radiative multipole moments, and these will be valid in weak-field, 
slow- motion situations. To formulate this approximation precisely, we introduce a characteristic radius r c , to be 
thought of as the radial coordinate of a typical portion of the source. [In other words, T; m (w;r) is assumed to be 
appreciably different from zero only for values of r comparable to r c ] We introduce also a characteristic time 1/uj c , 
to be thought of as the typical time scale over which the source moves. [In other words, T; m (w;r) is assumed to be 
appreciably different from zero only for values of u> comparable to lo c .] Finally, we introduce a characteristic velocity 
v <C 1, which will be the smallness parameter of our approximation. In terms of these quantities, the requirement 
that the source motions must be slow translates to 

io c r c = O(v). (3.13) 
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The virial theorem for gravitationally bound systems then implies that the gravitational field must be weak inside 
the source: 



M/r c = 0(v 2 ). 



(3.14) 



Finally, the slow-motion approximation implies that the scalar waves produced by the source's motion must have low 
frequencies: 



Mlo c = 0(v 3 ). 



(3.15) 



Our calculation of the radiative multipole mo men ts will be carried out t o or der v 3 beyond the leading-order expressions. 

We now proceed. It is evident that Eq. ( |3.6[ ) is nothing but Eq. (2.1), the generalized Regge- Wheeler equation, 
with s = 0. So we have, immediately, 



Equation ( 2.19 ) may therefore be substituted into Eq. ( 3. 10| ) . After using Eq. (3.4), we obtain 

4tt 



(3.16) 



Zlm{<Jj) 



^ ^ e 2iMa ) (ln4M| W |- ft )( 1 + 7rM | Ci; |)(_ ici;) i 



2(21 + 3) r + 1 [ 



p(u>, x) r l Y lr , 



i dx. 



where 



A = ^(l + l) + 



(3.17) 



(3.18) 



and dx = r 2 drdfl. 

The physical interpretation of this result comes more easily if we first invert the Fourier transform. The w-dependent 
prefactors complicate this procedure slightly, but an explicit expression for Zi m (u) can nevertheless be found (see 
Appendix D). We obtain 



Zi m {u) = Z lm (u) + 2M 



, . u — u . 



where a dot indicates differentiation with respect to u' , and 



Zi m {u) 



-Itt 



(2Z + 1)!! 



[rdu) 
2(21 + 



— -I — 

3) r 



0(v 4 ) 



Z[ m (u) du', 



p(u,x)r l Yi m (6,4>)dx. 



(3.19) 



(3.20) 



Equation ( 3.20 ) indicates that to leading order, Zi m (u) is given by the Z-th ret arded -time derivative of J pr l Yi m dx. 
This justifies our referring to these quantities as multipole moments. Equation ( 3.20] ) also shows that all corrections 



of order v appear inside the spatial integral, and that they depend on the detailed behavior of the source function. 
These corrections are near-zone corrections, and they are purely local in time: as with the leading-order term, they 
involve the value of the source at the retarded time u only. This is not so for the corrections of order v 3 , as is shown 
by Eq. (3.19): these involve the value of the source at all times prior to the retarded time u. Furthermore, the 



0(v ) corrections appear outside the spatial integral, and they are independent of the detailed behavior of the source. 



These are wave-propagation corrections. Equations ( 3.19| ) and (3.20), with /3i given by Eq. (3. IS), are equivalent to 
Eqs. 

We recognize the important distinction between near-zone and wave-propagation corrections. Near-zone corrections 
depend on the detailed behavior of the source and are local in time. Wave-propagation corrections, on the other hand, 
do not depend on the detailed behavior of the source, and are nonlocal in time. This nonlocality is heuristically 
understood as arising from the scattering of the radiation by the spacetime curvature surrounding the source. This 
scattering causes part of the information about the state of the source to be delayed further than what is strictly 
required by causality. The integral term in Eq. (3.19) is often called the tail term, and wave-propagation corrections 
are often called tail corrections. 
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IV. ELECTROMAGNETIC RADIATION 



A. Teukolsky equation 

In this section, we derive expressions for the radiative multipole moments of an electromagnetic field in Schwarzschild 



spacetime. The Maxwell equations for this spacetime have been cast, by Teukolsky 34|, in a form convenient for our 
purposes. In the Teukolsky formalism, the radiative part of the electromagnetic field is represented by the complex 
quantity $2 = F a/3 fh a np, where F a P is the field tensor, n a — —^(Z, 1,0,0) a null vector pointing radially inward, 
and fh a — (0, 0, r, — ir sin0)/v2 a spatial vector with zero norm. As before, a bar denotes complex conjugation. 
The field $2(x) has spin- weight s = — 1 (see Appendix A), and it can be decomposed according to 

$ 2 (x) = ^ J duY / Ri™(u;r)^ 1 Y lm (8,cf > )e-^ t , (4.1) 

Ira 

where -iYi m (9, 4>) are spherical harmonics of spin-weight —1 (see Appendix A). The sums over I and m are restricted 
by I > 1, \m\ < I. The radial function then satisfies the inhomogeneous Teukolsky equation |34j| , 



2f d * 1 



[wrf - 2icur(l - 3M/r) - 1(1 + 1) >Ri m (u;r) = T lm {u;r), (4.2) 



where / = 1 — 2M/r. 

The source term, T; m (w;r), is constructed as follows from J a (x), the (unspecified) current density. One first forms 
the contractions 

a J = -J a n cn _ X J = -J a fh a , (4.3) 

then evaluates the Fourier transforms 

s J(u, *) = ^J sJ(t, x) e 1 ^ dt, (4.4) 

and takes the projections 



s J lm (uj\r)= J s J(LO,x) s Y lm (6,4>)dQ. (4.5) 

The source term is finally given by p3 

T lm (oj;r) = 2ir^ s pi S T> s Ji m (uj;r), (4.6) 

s 

where the sum runs from s — to s = — 1, 



and 



with C = fd/dr + iuj. 



fvWTi) , = o_ i) (47) 



B. Solution 

Equation ([4.2[) is integrated by means of a Green's function, in a manner similar to what was done in Sec. III. We 



introduce two functions, Rf(ui;r) and Rf°(u>; r), solutions to the homogeneous Teukolsky equation [Eq. (4.2) with 
Ti rn (uj;r) = 0], with asymptotic behavior Q 
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R?(u;r^2M)~f e - iur ", 
R?(u,;r^ oo) ~ Qj n (o;)(iwr)- 1 e- iwr * + Q? ut (cj)(iur)e iur ' , 
i?;°°(w;r -> oo) - (W)e lwr *, 
where r* = r + 2M\n(r/2M — 1). In terms of these, the radial function is given by 

1 



(4.9) 



2iuQ l { l (uj) 



2M 



r /2 J/ 

where /' = 1 - 2M/r'. 

As before, we are mostly interested in the behavior of the radial function near r — oo. Equation ( 4.1 0| ) gives 



where 



R lm (cj; r ^ oo) ~ Z(cj)(itjr)e iur , 

T lm (w;r)R? (uj;r) 



dr. 



2iuQ?(u) J 2M r 2 / 

These quantities are the (Fourier transform of the) radiative multipole moments of the electromagnetic field. 



(4.10) 

(4.11) 
(4.12) 



C. Adjoint operators and Chandrasekhar transformation 



Equation ( 4.12j ) will not be our final expression for the radiative multipole moments. In Eq. ( 4.12j ), Ti m (uj\r) is 
obtained by applying the differential operators S T> on s Ji m (ui; r), the projections of the current density. To have to 
take derivatives of these functions is an inconvenience, and we would like to express the moments directly in terms of 
s J; m (w; r). It is easy to show, by straightforward integration by parts, that if we define the adjoint operators 



5 £r -l 



s = 
s = -1 



where C = fd/dr — iu>, then Eq. (4.12) is equivalent to 

(Ld) J 2 M 



s Ji m (u);r) s T>i Rf 1 (w;r) 
r 2 f 



dr. 



(4.13) 



(4.14) 



Although Zi m (uj) is now expressed directly in terms of s Ji m (Lo; r), Eq. (4.14) still will not be our final expression for 
the radiative multipole moments. We now want to write Rf 1 (u>; r) in terms of X, (w; r), the solution to the generalized 
Regge- Wheeler equation (with s = 1) considered in Sec. II. The relationship between these functions was trivial in the 
case of scalar radiation [cf. Eq. (3.16)]. That such a relationship exists in the case of gravitational radiation was shown 
by Chandrasekhar |35| , who also provided it explicitly. We show in Appendix E that for the case of electromagnetic 
radiation, the Chandrasekhar transformation is given by 



Q\ a (uj) 1(1 + 1)'" -4j n M 



rL ■ 



(4.15) 



Substituting this into Eq. ( 4.14 ), and taking into account the fact that CC = 1(1 + 1)/ /r 2 when acting on X l -(ui; r), 
we arrive at 



Zl m (uj) 



-2n 



poo 

^sPi / s Ji m (^;r) s TX"(cj;r)dr, 
JIM 



(4.16) 



1(1 + l)iLuA\ n (uj) 
where we have introduced the operators 

r 2 /- 1 /; 

-l(l + l)r 

We emphas ize th at A T, g (q;;r) den otes the solution to the s = 1 generalized Regge- Wheeler equation with boundary 
conditions (2.3)— (2.4). Equation (4.16) will be our final expression for the radiative multipole moments. 



s = 
s = -1 



(4.17) 
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D. Vector potential 



The physical meaning of the quantities Zi m (u>) becomes more transparent if we use Eqs. (4.1) and (4.11) to construct 
A™ d (x) , the vector potential describing the radiative part of the field. [This is defined similarly to $ ra d (x) , encountered 
in Sec. III.] 

We may choose a gauge in which A™ d (x) is purely transverse to the direction of propagation, which, at large 
distances from the source, is radially outward. This implies that the vector potential may be expressed as 



A™ d = At 



Ai 



(4.18) 



where m a = (0, 0, r, ir sin 9)/y2 is complex conjugate to m a . It should be noted that the quantity A(x) is complex, 
but that A™ d (x) is real. 

The radiative part of the electromagnetic field tensor, F"J^(x), can easily be computed from Eq. (4.18), and the 
Teukolsky field, $2(f) = ^rad 7 ™" 7 ^' follows immediately. Keeping in mind that we are working near r = oo, we find 



dA 

$2 = -A, a n a = - — , 
ou 



where u — t — r* . Combining this with Eqs. (hi) and ( 4.11 ) yields 

A(t,x) = - 53z Im (u)_iH m (d,i 



hn 



where 



Zi m {u) 



Z lm (u)e^ u dio. 



(4.19) 



(4.20) 



(4.21) 



This shows that Zi m (u) are indeed the radiative multi polc m oments of the electromagnetic field. 

The vector potential is obtained by substituting Eq.(4.2C) into Eq. ( 1.18] ). The spin- weighted spherical harmonics 
then combine w ith th e vectors m a and fh a to form the vectorial spherical harmonics described in Appendix A. [See 
Eq. (A16); Eq. JA15[) must also be used.] We find 



A^{t,x) = - r Y^ lm {u)Y^ 

1 771 



where we have introduced the charge multipole moments, 

1 



l lm (u) = — Z lm (u) + (-l) m Zi t - m (u) 
V2 L 



and the current multipole moments, 



i 

71 



Z lm (u) - (-l) m Zi,- m (u) 



(4.22) 



(4.23) 



(4.24) 



The reason for u sing this terminology will become clear below. For the time being we may mention that the first group 
of terms in Eq. ( 4.22| ), that involving the charge moments, has electric- type parity, while the second group, involving 
the current moments, has magnetic-type parity (see Appendix A). The fact that two sets of multipole moments are 
needed to form A™ d (x) is related to the fact that the electromagnetic field possesses two radiative degrees of freedom. 



E. Slow-motion approximation 



We now compute the radiative multipole moments ( 4.16|) in the slow motion approximation. The calculation is 
similar to the one presented in S ec. Il l C. 

We begin by substituting Eq. ( 2.19 ), with s 
obtain the lengthy expression 



1, into Eq. flll6D . After using Eqs. @, (gj), (|J), and ( |4.17| ), we 
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Zl m {u)) 



4tt 



(2Z + 1)!! 

itor 



l + l 
21 



l + l 
1 



%(lo) (-iu) 1 

,__M! 

2(21 + 3) Z 
(wr) 2 Z 2 - 1 M 



1 - 



____ M -(Z-l)- 



Z 2 - 21 - 1 M 



2(21 + 3) 



I 



0(v 4 



+ 0(v 4 )j (fp+J r ) r l % 
_i J r' -lYim dx 



dx 



Here, p(uj,x) = J*(o->,a;) is the Fourier transform of the charge density, and 



where /3j is given by Eq. ( 2.18 ) with s = 1: 



A = ^(Z + l) + 



l 



(4.25) 



(4.26) 



(4.27) 



2 2Z(Z + 1) 

Inspecting Eq. ( 4.25| ), we notice that it does not have the same mathematical structure as Eq. ( 3.17 ), which gives 
the radiative multipole moments of a scalar field. In particular, we see that Zi m (w) possesses correction terms that are 
linear in v [the terms iu>rp/(l + 1) and J r , the latter being one power of v smaller than p], as well as many third-order 
terms that depend explicitly on r, and which cannot be taken outside the integral. This apparently contradicts our 
expectation that Zi m (uj) should come with only near-zone corrections of order v 2 , and wave- propagation corrections 
of order v 3 . 



However, expression ( 4.25 ) is not unique, and we may use the continuity equation J a . a — to remove the unwanted 
terms. When written out explicitly, this reads 



4( r2jr L--i:(^ J + M, 



r* - ' ' T y/2r 

where d and d are the "edth" differential operators described in Appendix A, and 

i J = —J a m a . 



(4.28) 



(4.29) 



The continuity equation gives rise to an integral identity if we multiply both sides by r™ aYi m (9, (f>) and integrate over 
dx. After a Fourier transform and several partial integrations [using Eqs. (A13) and (A14)], we obtain 



iw I pr n Yi m dx = n I J r r n 1 Yi m dx 



1(1 + 1) 



_1 J -lYi m - 1 J lYi m r 



dx. 



(4.30) 



We now use this identity to remove all terms proportional to iuip(uj,x) in Eq. (4.25). After some remarkable cancela- 
tions, wherein all unwanted terms disappear, we arrive at our final expression for the radiative multipole moments: 



2>im(w) 



4tt 



(21 + 1)!! 

+ 0(v 4 



Tt(u) (-iuj) 1 



l + l 
21 



1 - 2 ( z + i)(2z + 3) M 2 -a-D- 



J r ) r l a Yi m dx- - 



0(v 4 



luir 
T+l 



-\J-xYi m + \ J\Y lm r dx 



1 - 



(tor) 2 



I 2 -I M 



2(2Z + 3) Z 



(4.31) 



We see that this expression has the expected form, with all 0(v 2 ) corrections occurring inside the spatial integrals, 
and all 0(v 3 ) corrections occurring outside. 

We now separate Zi m (u) into charge and current moments, according to the Fourier transform of Eqs. ( 4.23 ) and 
(Ol. This gives 



+71 



(21 + 1)!! V Z 

-(Z-l)- + 0(« 4 ) 
r 



Ti(u>) (-iuj) 



1 + 3 



2(Z + l)(2Z + 3) 



(cor) 2 



_ iLur - 



r l Y lri 



i dx 



(4.32) 
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and 



(2Z + 1)!! 



+ 0(v 4 ) 



(urf 



2(21 + 3) 



I 



1 M 

r 



-iJ(u,x)^ 1 Yi m (9,(f>) + iJ(w,5)iY/ m (0,< 



r dx. 



The corresponding expressions in the time domain (see Appendix D) are 



ll m (u)=I lm (u)+2M 



and 













+ A +7 


/ — oo 


V 4M / 



5 im (u) = 5 iro (tt) + 2M 

where a dot indicates differentiation with respect to u' . We have defined 

I 



Ii m (u') du' 



Si m {u') du', 



Ilm(u) 



47T 



l + l I d_ 

da 



(2Z + 1)!! V I 



r 



1 



1 + 3 



2(l + l)(2l + 3) 



(rd u ) 2 



p(u, x) J T (u, x) 



r l Y lm {e,<l>)d£ 



and 



2V2i 



(21 + 1)!! \du 
iJ(u,x) -ifi 



(rd u 



I 2 - 1 M 



2(21 + 3) I 
+ iJ(u,x) iYi m (9 



0(v 4 



r l dx. 



(4.33) 



(4.34) 



(4.35) 



(4.36) 



(4.37) 



Equations ( 4.34 )-( 4.37 ) have the same mathematical structure as Eqs. ( 3.19 ) and ( 3.2C ), which give the radiative 
multipole moments of a scalar field. The physical meaning of these equations is therefore exac tly th e same as in 
Sec. Ill, and the discussion appearing at the end of Sec. Ill E need not be repeated. Equations (4.34)-(4.37), with 
[3i given by Eq. (4.27), are equivalent to Eqs . (|l.8[ )-( 1.12| ), once the spin-weighted spherical harmonics have been 
converted into the vectorial harmonics of Eq. (|A16|). 



V. GRAVITATIONAL RADIATION 



A. Teukolsky equation 

In this section, we derive expressions for the radiative multipole moments of a gravitational-wave field in 
Schwarzschild spacctime. Specifically, we consider a tensor field h a p(x), defined as the difference between the metric 
of the perturbed spacetime and the Schwarzschild metric. Field equations for h a p(x) are obtained by linearizing the 
Einstein equations for the full metric. It is therefore assumed that the perturbation is small. Teukolsky |34| has cast 
the field equations for h a p(x) in a form convenient for our purposes. We briefly summarize this formulation here. 

In the Teukolsky formalism, the radiative part of h a p(x) is represented by the complex- valued function ^4 = 
—C a /3 1 sn a fh^n 1 fh s , where C a j3 7 s is the perturbed Weyl tensor, and n a and fh a are defined as in Sec. IV. The field 
^4(0;) has spin-weight s = —2 (see Appendix A), and it can be decomposed according to 

* 4 = ^ J diuJ2Ri m (^r)~2Yi m (8,<j>)e- luJ \ (5.1) 

lm 

where -2^1™ ($, 4 1 ) are spin-weighted spherical harmonics (see Appendix A). The sums over I and m are restricted by 
I > 2 and |m| < I. The radial function then satisfies the inhomogeneous Teukolsky equation JmJ, 
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dr 



Rim{u;r) = T iTO (w;r), 



(5.2) 



where / = 1- 2M/r and U{u)\r) = /^[(wr) 2 - 4iw(r - 3M)] - (I - l)(/ + 2). 

The source term to the right-hand side of Eq. ( |5.2| ) is constructed as follows from T a/3 (x), the (unspecified) stress- 
energy tensor responsible for the perturbation. The first step is to form the contractions 



One then evaluates the Fourier transforms 



_ 2 T = T a pm a m 



and takes the projections 



Finally, T lm (u; r) is given by p 



s f(uj,x) 



1 

2^ 



s T(t, £) e iaj * dt 



s fi m (w;r) = / s T(u,i) a Yi 



T lm (uj;r) = 2ir*^2 s pi S V s fi rn (uj; r), 



(5.3) 



(5.4) 



(5.5) 



(5.6) 



where 



2^/ (2 + !)( ; + 2) s = 
.Pi = <( 2 v /2(/ - 1)(Z + 2) s = -1 

1 s = -2 



(5.7) 



and 



s = 



rfCr^f^Cr 



(5.8) 



Here, C = fd/dr + iui. 



B. Solution 



The inhomogeneous Teukolsky equation ( |5.2p can be integrated by means of a Green's function, in a manner 
similar to what was done in Sec. IV. Here also, the form of the radial function can be simplified by introducing adjoint 
operators s Z)t, an d by expressing it in terms of X^(lo\ r)/^4J n (w), wher e X,( uj; r) is the solution to the Regge- Wheeler 
equation — Eq. (2.1) with s = 2 — with boundary conditions (2.3), (2.4). These manipulations are described in 
detail in Ref. plJI , and they will not be displayed here. The conclusion is that at large distances, the radial function 
is given by 



2 -5- 



(lo) re 



3 „iojr* 



(5.9) 



where 



Am 



-2tt 



iwKi(w)^ n (w) 



sPl 



rf- 2 s f lm {uj;r) S T X?(w;r). 



(5.10) 



2AI 



Here. 



Ki{u) = ~ (I - 1)1(1 + + 2) - 12iMuj 



(5.11) 



and 
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r = 2(1 - 3M/r + iur)rf — + f\l(l + 1) - 6M/r) + 2itor(l - 3M/r + iur), 
-iT = -/{ [1(1 + 1) + 2iur]rf^ + 1(1 + l)(f + iwr) - 2(^r) 2 }, 
_ 2 r = / 2 {2[(Z - + 2) + 6M/r]rf-^- + (I - l)(l + 2) [1(1 + 1) + 2iujr] + 12/M/r}. 
The quantities Zi m (uj) are the multipole moments of the radiative part of h a p(x). 



(5.12) 



C. Metric perturbation 

The gravitational- wave field, h r ^(x), can be obtained from the behavior of ^(x) at large distances Choosing 
the 9 and 4> directions as polarization axes, the two fundamental polarizations of the gravitational waves are given by 

h + (t, x) - ih x (t, x) = - y~] Z lm (u) -iY lm (6, 4>), (5.13) 



l in 



where u = t 



and 



Zim(u) 



(5.14) 



This shows that Zi m (u) are indeed the multipole moments of the radiative field. It should be noted that while these 
quantities are complex, h + (x) and h x (x) are real. 

In the traceless-transverse gauge, the gravitational-wave tensor is given by 



= (h+ - ih x )m a mp + (h+ + ih x )m a m l 3, 



or, after substituting Eq. (5.13), 



(5.15) 



(5.16) 



Here, h r ^(x) is expressed in terms of the tensorial spherical harmonics described in Appendix A. The mass multipole 
mom ents T ; m (u) , and the current multipole moments Si m (u) , are related to Z\ m (u) by the same equations as Eqs. ( 4.23| ) 
and (fl24|). 



D. Slow-motion approximation 

We now calculate the radiative multipole moments in the slow-motion approximation. We proceed as in Sec. IV E. 
Substituting Eq. (§1J|), with s = 2, into Eq. ( pU0| ), and using Eqs. (|J), @, ( P^ and (|I|), yields 



Zim(u)) 



16tt 



— Ti(u) (-iu;) 1 J2sPim(u) 



where 



with 



(21 + 1)! 

Ti(u) = e 2lM ^ n4M ^-^(l + irM\Lu\) 



Mi = Pi 



(l-l)l(l + l)(l + 2Y 



while Pi is given by Eq. (2.18) with s — 2: 



a = + 1) + 1 - 1 



2 1(1+1) 



(5.17) 



(5.18) 



(5.19) 



(5.20) 
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We have also introduced 



o-P;m(w) — 2-t 



/ 0+1)0 + 2) 
' (i-l)' 
Z + 4 
2(Z + 2)(2Z + 3) 



1 - 



1 + 9 



L-Pzm(w) = -2 



2(Z + 2) 
I - 1 



1 - 



2(Z + 1)(2Z + 3) 

Z 3 + 3/ 2 + I - 4 M 
Z(Z + 2) r~ 

Z 2 + 3Z + 6 



(wr)^-(Z + 2)— + 



M 2iujr 



l+l 



2Z(Z + 1)(2Z + 3) 



(wry 



+ O(v^^r l Y lm dx, 

(I- 1)0 + 2) M iwr 
F + ~T 



I 2 + 31 + 6 



2(Z + 1)(Z + 2)(2Z + 3) 



K) 2 " 



I 3 + 3l 2 



M 



+ 1)0 + 2) r 



+ 0(v A ) 



iT 



-2-PiroM = 



Z + 2 



i 



+ 1)0 + 2) 



+ 0(?; 2 ) 



_ 2 iV - 2 Yi m dx. 



(5.21) 

(5.22) 
(5.23) 



To better keep track of the relative importance of each term in Eq. (5.17), we decompose the stress-energy tensor 
according to 



03 = T tr , 
oP = T rr , 
t = T af3 m a mi3, 



P = T tt , 

-13 = -T ta m a , ij = -T ta m a , 

-lP = -T ra fh a , ip = -T ra m a , 
-2t = T a ^fh a mj3, it = T a ^m a mp. 



(5.24) 
(5.25) 
(5.26) 
(5.27) 



Thus, if p is cons idered to be a quantity of order unity, then s j = 0(v), s p = 0(v 2 ), and s t = 0(v 2 ). Inspection of 
Eqs. (pll)-(p3D - in which we substitute Q T/f 2 = \{p + 2 Q j / f + oP / J 2 ), _ x T/f = |(_ij + _ 1 p//) ) and _ 2 T = _ 2 < 
— then reveals that as it stands, Zi m (uj) contains many unwanted terms of the sort encountered in Sec. IV E: terms 
which are first order in v, and 0(v 3 ) terms that cannot be pulled outside of the spatial integrals. 

In Sec. IV E, the unwanted terms were removed by invoking the continuity equation, J a a — 0. Here, they are 

removed with the help of the energy-momentum conservation equations, T"^ = 0. When written out explicitly, 
these become 



o3, t 



r 2 \Vf 



op 



+ -L- (d-ip + d lP ) + hfp - y~ t = o, 

V2r '2 r 



-d,t + ^(r 3 -iP), r + ^ b - 2t + B ot) = 
d,t + ^{r\p) r + -±-(£> t + d- 2 t) = 0, 



(5.28) 
(5.29) 
(5.30) 
(5.31) 



where d and d are t he "ed th" differential operators described in Appendix A, and /' = df jdr = 2M/r 2 . 

Equations (|Bl^-([5~3i"l) give rise to a number of integral iden tities , whic h we write in the frequency domain, and 
which are easily established by partial integration, using Eqs. (A13) and (A14). We shall need the following two 
identities, which involve no approximation: 



iuj I pr n Yi m dx = 



{n-rf/f) ~ 3 r n - l Q Y lm dx 
-13-1% 



1(1 + 1) 



ijlYlm) 



dx 



(5.32) 



and 



-ijr n -iYi m dx = -(n - I) 



-ipr n 1 -iYi m dx 



0-1)0 + 2) 



- 2 Yi m dx 



1(1 + 1) 



otr n 1 Y lm dx. 



(5.33) 
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We shall also need the following identities, which are valid in the slow-motion approximation: 

(_iJ-iYi m - iJiYim) r l dx 











j iivr 




Mi- 


hi) 



T (_ip -\Ylrn - lP iYi m ) r l dx, 
where terms of order v 5 and higher (with p taken to be of order unity) have been discarded, and 
[iMtu + 0(v 5 )] -!jr l -lYhndaS =-(l-l) [ [M/r + 0(v 4 )] -ipr 1 ^Y lm dx, 



(5.34) 



(5.35) 



which follows from Eq. ( |5.33 ) after multiplying both sides by M. 



These identities are used to remove all unwanted terms from Zi m (uj), as given by Eqs. (5.17)— (5.23). After a rather 
long calculation (which spans several pages), we arrive at the following expression for the radiative multipole moments: 

8tt 



(21 + 1) 



■Ti(u>) (-iu) 



I J / (f + l)(f + 2) 
(I - 1)1 



1 



1 + 9 



2(Z + 1)(2/ + 3) 



(wry 



(l + 2)— + 0(v 4 ) 
r 



p + Q p + 2 a t + j-^-j qj ) r Yi m dx 



2(1 + 2) 
l-l 



1 + 4 



2(l + 2)(2l + 3) 



( „. r)2 (l-m + 2)M 



24iMu 



0(v 4 



-i3- 



-ip r -iY lm dx 



r (l-l)l(l + l)(l + 2) 
I 

1 



2(1 + 2) 
l-l 



2(l + 2)(2l + 3) 



(corf 



(Z-l)(Z + 2) M 

l 7 



+ 0(v 4 ) 



icur 

lj+ lP ) r iYi m dx 



(5.36) 



Notice that the second and third integrals differ by a term proportional to iMuj. We see that this expression for 
the radiative multipole moments has the expected form, with all 0(v 2 ) correction terms occurring inside the spatial 
integrals, and all 0(v 3 ) corrections occurring outside. 

We now separate Zi m (u>) into mass and current moments, according to the Fourier transform of Eqs. ( 4.23| ) and 
(4.24). (As was pointed out in Sec. IV C, those equations are valid also in the case of gravitational radiation.) We 
find that the mass moments are given by 

16vr 1(1 + l)(l + 2) 



Ilm(w) 



(21 + 1)!! V 2(1-1)1 
+ 0(v 4 ) 



Ti{u) (-iu) 1 



2(l + l)(2l + 3) ^ - {l + 2) T 



Aiuor 



p(uj, x) + p(u, x) + 2 t(u), x) + j(u, x) 



r oYi m (6,4>) dx, 



(5.37) 



and that the current moments are given by 

~ lQiri 
Sl m (w) 



0(v 4 



(21 + 1 

(Z-l)(Z + 2) M 

I 7 



1 + 4 



2(l + 2)(2l + 3) 



(wrf 



-1%, 



_. nor _. „ 

ij(u),x) + -—— lP(U),X) 



1 + 2 



{Ylr. 



r dx, 



(5.38) 



where 



T l i (w) = e 2lM ^ lniM ^- l3 ''>(l + TTM\u;\). 



(5.39) 
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Notice that different cons tants (/ij fo r the mass moments, Pi for the current moments) appear in %{ui) and 7^ (w); 
these are defined by Eqs. and fTl]). 

The corresponding expressions in the time domain (see Appendix D) are 



Tlmiu) = Ii m (u) + 2M 



and 



In 



In 



AM 



AM 



+ (J>i + 7 



+ Pi +7 



S im (w) = Sim(«) + 2M 
where a dot indicates differentiation with respect to v! . We have defined 

Z + 9 



hm{u') du 



Slm{u') du' 



16tt 
(21 + 1)11 

+ <3(t> 4 ) 



' (Z + l)(Z + 2) 
2(Z — 1)Z Vdu 



1 



2(Z + 1)(2Z + 3) 



(rd u ) 



p(it, x) + op(u, x) + 2 t(u, x) 



Ard u 
l + l 



oj(u,x) 



I da? 



and 



Si m (u) = -- 



16ni 



1 + 2 



(2Z + 1)!! V I -I \du 

l T 



l + A 



2(/ + 2)(2/ + 3) 



{rd u f 



ij(u,x) 



ij(u,x) 



1 + 2 



ip(u,x) 



lYlr 



rd u 
1 + 2 



dx. 



-%p(u,x) 



-iY h 



(5.40) 



(5.41) 



(5.42) 



(5.43) 



The interpretation of these results is exactly the same as for the cases (scalar and electromagnetic radiation) considered 
previously. Equations (5.40)-(5.43), with (ii = ps rav ' cunen anc j ^ = ^grav.mass gj ven gq S ^ (5 20) and (5.19), 
respectively, are equivalent to Eqs. (1.14)-( 1.19| ), once the spin-weighted spherical harmonics have been converted 
into the vectorial harmonics of Eq. ( A16| ). 



VI. PHYSICAL ORIGIN OF THE TAIL TERM 



A survey of Sees. Ill, IV, and V reveals that for scalar, electromagnetic, and grav it ationa l radia tio n, th e tail 
correc tion to the radiative multipole moments takes the universal form [cf. Eqs. ( 3.19 ), ( 4.34 ), ( 1.35 ), ( 5.40 ), and 

(PD] 



Mi m {u) =M lm (u) + 2M 











U u - U ') 


+ ci + 7 


t —00 


[ \ AM J 



Mimiu) du. 



(6.1) 



Here, A4i m {u) stands for Zi m (u) in the case of scalar radiation, and for Ti m (u) and Si m {u) in the case of electromagnetic 
and gravitatio nal ra diation. Similarly, Mi m (u) stands for either Zi m {u), Ii m (u), or Si m (u). The constant q stands 
for Pi [cf. Eq. (|2.18| )], except for the mass multipole moments of the gravitational- wave field, for which c/ stands for 
m [cf. Eq. Q5.19DI. 



The physical interpretation of Eq. (f 

that while the correction terms of order v z that appear in Mi 
detailed behavior of the source, the correction terms of order v 3 ■ 



is clear, and was first given at the end of Sec. Ill C. Equation (6.1) shows 
2 that appear in 



(it) are near-zone corrections that depend on the 
- the terms under the integral sign, or tail terms — 
are due to wave-propagation effects, and are independent of the detailed behavior of the source. And while the 0{v 2 ) 
corrections are local in time, the 0(v 3 ) corrections introduce a nonlocality in the radiative multipole moments. This 
nonlocality is understood as ari sing from the scattering of the radiation field off the spacetime curvature generated 
by the mass M, and as Eq. (3.1) shows, the tail term is indeed proportional to M. 
The mass parameter appears in two places in the Schwarzschild metric, 



ds 2 = -fdt 2 + f~ l dr 2 



Hde 2 



sin 



2 ). 



(6.2) 
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where / = 1 — 2M/r. It enters in g u and in g rr , which are both involved in the calculation of the tail correction. 
Because the mass parameter is the same in both components of the metric, it is impossible to tell, on the basis of our 
previous calculations, whether it is gu that is "mostly responsible" for the tail effect, or whether it is g rr , or whether 
both components are "equally responsible" . In other words, our previous calculations cannot tell us how the temporal 
and spatial curvatures separately contribute to the tail effect. This is the question we now wish to examine. We shall 
answer it by artificially introducing an additional mass parameter in the description of our spacetime. For simplicity, 
we shall restrict attention to the case of scalar radiation. 

We consider a static, spherically symmetric spacetime with a line element of the most general form, 



ds 2 



-fdt 2 + g^dr 2 



2 (d0 2 



where / and g are two arbitrary functions of r obeying the following restrictions, 
asymptotically flat, so that the metric functions must behave as 

/ ~ 1 - 2M/r, g~l- 2 7 M/r, 



(6.3) 

First, the spacetime must be 



(6.4) 



at large distances (r ^> M). Here, M is the gravitational mass of the system responsible for the gravitational field, 
and 7 is a parameter that measures the failure of the metric to match the Schwarzschild form at large distances; 
7M can be thought of as the system's inertial mass, and the Schwarzschild behavior is recovered by putting 7=1. 
(This parameter has the same meaning as 7 in the parameterized post-Newtonian formalism p6[ . We nevertheless 
use the notation 7 to distinguish this quantity from the Euler number 7.) Second, we assume for concreteness that 
the metric describes a black-hole spacetime, so that both f(r) and g(r) vanish at a common radius r$. Regularity of 
the spacetime at the event horizon further demands that the ratio / jg be finite and nonvanishing at r = r$. (Our 
conclusions are insensitive to this second set of assumptions.) Apart from these requirements, f(r) and g(r) will be 
left unspecified. 

We consider the scalar wave equation, Q$>( x) = —4Trp(x), in a spacetime with line element (6J3). After separation 
of the variables, according to Eqs. (3.2), (3.4), and ( |3.5|) , the radial function is found to satisfy 



dr* 



1(1 + 1) 



f r 



-Rl m (u;r) = fTi m (uj;r), 



(6.5) 



where dj dr* — \[Jg d/dr, and a prime indicates differentiation with respect to r. This equation is integrated by means 
of a Green's function, constructed from two linearly indepen dent solutions to the homogeneous equation. These are 
denoted Rf(uj;r) and i? ; °°(tj; r), and are defined as in Eq. (3.7), with r* = J dr/y/Jg, and with r = replacing 
r = 2M. The solution at large distances is then given by Eq. (3.9), with 



Zi m {uj) = 



1 



r JJgTi m (uj-r) Rf(cj;r) dr. 



Substituting this into Eq. (3.2) yields 



rati 



(*, x) = ~y] Zi m (u) Y hn ( 



(6.6) 



(6.7) 



where Zi m (u) is the inverse Fourier transform of Zi m (ui) 
is now given by 



and u = t — r* is retarded time. According to Eqs. (6.4), 



at large distances, where 



r* ~ r + 2aMm(r/2M) 



1(1 + 7). 



(6.8) 



(6.9) 



Equation (5 



agrees with the Schwarzschild definition when 7 = 0" = 1. 
We now wish to calculate Zi m (oj), the radiative multipole moments, in the slow-motion approximation. 



The first 

step is to integrate the homogeneous version of Eq. (6.5) in the low- frequency limit. The calculation proceeds as in 
Sec. II and Appendix C, and uses the approximations (3.4) for f(r) and g(r); these steps will not be duplicated here. 
Defining e = 2Mu> and z = lot, we eventually find 
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= 2(-i) l + 1 e iae{ln2e - f,l \\ + \ae)z\ \l ~ aeAi(z)]ji(z) + aeBi[z)ni{z) 



(I + 1)7 — 0- . . . lj + a . , . 2x 



2(2/ + 1) 



2(2/ + 1)' 



where Ai(z) and S;(z) are denned by Eqs. ( 2.14 ) and ( 2.15 ), respectively, and 



Pi = 1>(l + 1) + 



2a' 



(6.10) 



(6.11) 



Evaluating Eq. ( 6.1C ) for z <C 1 yields 

ffifor) _ 2 



2i CT M W (ln4M|a,|-ft) A + ^/U) (_ ia;r )l+l 

Q5 n (a;) (21 + 1)!! 1 
f (wr) 2 



x < 1 



(Z + 1)7-0- — + <3(t> 4 
r 



2(2/ + 3) 

The second step is to substitute Eq. ( |S.12 ) into Eq. (|6^), using Eq. (6.4) once more. We obtain 



(6.12) 



Zlm(u) 



An 



' e 2i*Mu{UiM\"\-Pi)( 1 + ir aM\u\)(-iu) 
(2Z+1)!! 1 

p(u>,x)r l Yi m (6,4>) dx 



(turf 



2(21 + 3) 



( 2 /-l)7+lM + ■ 
2 r 



The corresponding expression in the time domain is 
Zi m (u) = Z [m (u) + 2aM 



, / u — u 1 



Zi m (u') du', 



where 



-In 



(2Z + 1)!! 



1 



(rd u ) 



(21 -1)7+1 M 



2(2/ + 3) 



0(v 4 ) 



p(u,x)r l Yi m (9,cj>)dx. 



(6.13) 



(6.14) 



(6.15) 



These are t he ra diative multipole mom ents o f a scalar field in a spac etime with line e lement ( |6.3|) , ( |6.4| ). Equations 
and ( |6T5| ), with A given by Eq . ( |6.1l] ), are equivalent to Eqs. ( |l.20| ) and ( |l.2l| ). 
We see from Eqs. ( 6.14 ) and ( 6.15 ) that 7 and a both appear in the near-zone and wave-propagation correction 
terms. In particular, the tail integral is now proportional to a = 5(1+7). This allows us to conclude that in 
general relativity, temporal and spatial curvatures contribute equally to the tail correction. This result is striking, 
because the same conclusion is known to hold in two other situations: the deflection and time delay of light by the 
gravitational field of a massive body. Indeed, in a parameterized post-Newtonian calculation [|6| , the deflection angle 
and the time delay are both found to be proportional to |(1 + 7). (Here, 7 is the parameter that measures how 
much spatial curvature is produced by a unit rest mass; it is equal to unity in general relativity.) The statement 
that temporal and spatial curvatures contribute equally therefore applies to two very different physical situations. 
While the deflection and time delay of light arc both high-frequency, geometric-optics phenomena, the tail effect is 
very much a low-frequency, wave-like phenomenon, and the discovery of such a similarity in such different situations 
could not have been expected on physical grounds. However, this similarity is not entirely surprising on mathematical 
grounds: The factor of a that appears in front of the tail integral is essentially the same a that appears in the new 
definition of r* , Eq. flS.Ej ); since radial light rays propagate along curves of constant t — r* or t + r* , it is perhaps not 
surprising that a should also appear in expressions for the deflection angle and the time delay. 



VII. SPACETIME APPROACH 



A. DeWitt-Brehme Green's function 

The mathematical methods employed in the previous sections of this paper to derive expressions for the radiative 
multipole moments of integer-spin fields were based upon a separation of variables approach made possible by the 
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symmetries of the Schwarzschild solution. Although the physical interpretation of our results, in terms of near- 
zone and wave-propagation corrections, is quite clear, we cannot claim that the physical picture is particularly well 
represented by the mathematics involved in bringing the problem to a solution. Indeed, the physical meaning of our 
expressions became clear only after performing the inverse Fourier transform that gave the multipole moments in the 
time domain; by themselves, the frequency-domain expressions did not have a very compelling interpretation. 

In this section, we offer an alternative derivation of the radiative multipole moments in which the mathematics 
reflects the physics every step of the way. For simplicity we shall again restrict attention to the case of scalar radiation. 

Following the seminal work by Hadamard (!]], DeWitt and Brehme || considered the scalar wave equation, 0<f>(x) — 
—4irp(x), and its Green's function satisfying DG(x, x') — —S(x, x'), in an arbitrary spacetime with metric g a p. These 
equations imply that the scalar field can be expressed as 



>$>{x)=Att J G(x,x') p(x')dx', 



(7.1) 



where dx' = \g(x r )\ 1 ^ 2 d 4 x', with g — det(g a p). Assuming that the field point x belongs to the normal convex 
neighborhood of the source points x', DeWitt and Brehme found that the retarded Green's function takes the form 



G(x,x') = i- 



(x, X 



u(x, x') 8(a) — v(x, x') 6(—<j) 



(7.2) 



Here, a(x, x') is the world function first introduced by Synge p7| , and equal to one-half the squared geodesic distance 
between x and x'; a is positive if the points are spacelike related, negative if the relation is timelike, and zero if x 
and x' are joined by a null geode sic. The functions u(x,x') and v(x,x') are nonsingular in the limit a — > 0, and 
are obtained by substituting Eq. ( |7.2|) into the differential equation for the Green's function. Finally, 0(x,x') is a 
time-ordering function, equal to unity if x is in the causal future of x' , and zero otherwise. 

As can be seen from Eq. (7.2), the retarded Green's function splits naturally into a direct part (the first term), 
which has support on, and only on, the past light cone of x (all points x' such that a = 0), and a tail part (the second 
term), which has support inside the past light cone (all points x' such that a < 0). This, in turn, implies that $(:c) 
will also be split into direct and t ail p arts, as w as observed in Sec. III. We therefore see that contrary to our previous 
mathematical formulation, Eqs. (7.1) and (7^2) reflect the physical picture quite closely. 

In the remainder of this section, we calculate G(x,x') in a weak-field approximation (relying on previous work 
by DeWitt and DeWitt |p8[ ), and derive an expression for the radiative multipole moments of the scalar field. Not 
surprisingly, our answer will agree with what was obtained in Sec. Ill, Eqs. ( 3.1 9| ) and (3.20). Although this calculation 
tells us nothing new in terms of the final answer, it is still instructive, because of the fact that the mathematical origin 
of the tail correction is clear from the outset — it follows directly from the tail term in the Green's function. 



B. Direct term 



We begin with the calculation of the direct part of the field, 

$directO) = / 0(x, x 1 ) u(x, x') 5(a) p(x') dx' . (7.3) 



This involves the evaluation of o(x, x') and u(x, x'). We shall work in the weak-field approximation (that is, linearized 
gravity in harmonic coordinates), and express the metric as g a @(x) = r) a p + h a p(x), where r] a /3 is the metric of flat 
spacetime in Cartesian coordinates, and 

h a p(x) = (7-4) 
\x\ 

Here, 5 a p is the Kronecker delta, and for any three-vector s, |s| 2 = s ■ s = S a bS a s b . It is assumed that the source of 
the gravitational field is a point mass located at the origin of the coordinates. 
The world function is given by [ |39| 

. 1 f d£ a d^ „ ,„ c , 

C r(x,x) = -J c g a , — —dX, (7.5) 

where C is the geodesic relating the points x' and x, £ a (A) the equation o f th is geodesic, and A an afflne parameter on 
the geodesic, normalized so that £ Q (0) = x' a and £ Q (1) = x a . Equation (7.5) follows immediately from the geometric 
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meaning of the world function. Because Eq. (7.5) is an action principle for the geodesic equation, an error of order e in 
the specification of C is translated into an error of order e 2 in a(x, x'). Since we wish to evaluate a(x, a; 7 ) accurately 
to first order in the formally small parameter M, it is sufficient to approximate C by the straight path |3£j 



f* (A) = x' a + X(x a -x' a ). 



Substituting this into Eq. (7.5) and discarding all 0(M ) terms, we obtain 



1 



1 dX 



a(x,x') = --[ 1-2M / j)(f-t_)(f-i 



(7.6) 



(7.7) 



where 



f 1 d\ 

t±=t±\x-x'\±2M\x-x'\ -P 



and £ = |£|. Equation ( |7.7[ ) implies 



5(a) 



\x — x 



S(t' -t_) + 6(t' -t+) 



(7. 



(7.9) 



and the second term vanishes when 5(a) is multiplied by 9(x,x'). To evaluate the integral in Eq. (7\. 
Eq. (Q to write 

£ = a/^' 2 + 2Af • (z - £') + A 2 \x - x> | 2 , 



where r' — |x'|; we also define r = \x\ and n — xjr. The integration is elementary, and Eq. (7.£) becomes 

t ± =t±\x~ x'\±2M\n(2r/s'), 



where s' 



n ■ x' , and where terms of order unity have been discarded in the logarithm. 



we use 
(7.10) 

(7.11) 
(7.12) 

where R a p is the Ricci tenso r. B ecause R a p is already proportional to M, the geodesic C can once again be approxi- 
mated by the straight path ( |7.6| ). Now, R a p oc S(x ) in t he point-mass approximation, and there is only one path C 
which gives rise to a nonvanishing integral in Eq. (|7.12|) — the path for which x' and x are diametrically opposite. 
Because C forms a set of measure zero in the space of all paths connecting source points x' to a given field point x, 
the fact that u(x, x') ^ 1 for this path has no effect on $dircct(aO- Therefore, we can safely set 



In the weak-field approximation, u(x,x') is given by p9| | 

u(x, x') = l + \(x a - x' a )(x p - x' fJ ) I R a0 \(l - A) dX, 



u(x, x') = 1 



in the following. 



Substituting Eqs. (|7_9|), ( |7.1l| ), and ( |7.13| ) into Eq. (|7J) yields 



dx' 



where dx' = \g(x')\ 1 ^ 2 d 3 x = (1 + 2M/r') d 3 x. At large distances, this becomes 



4> 



rad 
direct 



where 



(t,x) = - I p[u + n- x' + 2M\n(s' /2c),x'] dx" , 



u = t-r- 2Mln(r/c) 



(7.13) 



(7.14) 



(7.15) 



(7.16) 



is retarded time, with c an arbitrary constant. This definition of retarded time is similar to the Schwarzschild 
expression, and c will eventually be chosen so that the two definitions agree. 
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We now invoke the slow-motion approximation and expand p in a Taylor series about u. (The approximation 
ensures that the series converges.) This gives 

1 00 1 r 

$^ r d cct (M) = -^- p^(u,x')[n-x' + 2MH S / /2c)] n dx', (7.17) 

r n=0 U ' J 

where p' n ' = d n p/du n . After discarding all terms of second and higher order in M, we arrive at 

1 00 . 

*d?r d ect(*, ^) = ~ £ 3 / ^ ^ f ') ^ ' df ' 

+ ^ - / ln(s'/2c)(n-x') n dx'. (7.18) 

r n=0 

This is our final expression for the direct part of the radiative field. 



C. Tail term 



The tail part of the scalar field is 

$taii(aO = - J 0(x, x 1 ) v{x, x') 6(-a) p(x') dx', 



(7.19) 



and it is now our task to evaluate this. 

An expression for v(x,x'), accurate to first order in M in a weak-field approximation, was derived by DeWitt and 
DeWitt pi, who find 



2M d 2 



vix, X 



\x-x'\ dt'dt 



9{r + r' +t' -t) In 



r + r + \x — x 



r + r — \x — x' 



+ 6{t -t 1 -r' -r) In 



t-tf + \x-x'\ 
t-t' - \x-x'\ 



For large r, this reduces to 



where s' = r' + ft ■ x! . 



v(x, x') = — 



2M 



5{u - u' - s') 0(u - u' - a') 
s' (u — u') 2 



u' = t' - r + \x - x'\ ~ t' - ft ■ x" ', 



(7.20) 



(7.21) 



(7.22) 



and u = t ~ r is retarded time. [The true retarded time is given by Eq. (7.16) and differs from t — r by a term 
2Mln(r/c). Nevertheless, u — t — r is the appropriate expression to use in the calculation of the tail term when 
working to first order in M, be cause v(x , x') is already proportional to M.j 

We now substitute Eq. (|7.21|) i nto Eq. Q7.19D , taking into account that 6(x, x') 6 (-a) = 9(t— \x-x'\-t') = 6(u-u'), 
according to Eqs. ( |7.7[ ) and (7^8). We obtain 



*28(t,*) = — 

r 



5{u-u'-s') e(u-u'-s') 



p(v! + ft ■ x' , x') du' dx' , 



(7.23) 



s' (u — u') 2 

where dx' — \g(x')\ 1 ^ 2 d 3 x' = [1 + O(M)] d 3 x' . After two partial integrations and a few lines of algebra, this becomes 

2M 



— p(u — r' , x) In s — I p\u — (,x) ln(C + ft ■ x) dQ 



+ / p{u + ft ■ x, 5?) ln(w — u ) du 1 



dx. 



(7.24) 
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where a dot indicates differentiation with resp ect to either u or 

The C-integral to the right-hand side of Eq. ( |7.24 ) can be evaluated explicitly if p(u — £, x') is expanded in a Taylor 
series about £ = 0. (Again, the slow-motion approximation ensures that this series converges.) This results in an 
infinite sum of terms involving the integrals J £™ ln(£ + n ■ x') dC t which can be expressed in closed form (Ref. [^o| , 
p. 205). After a rather long but straightforward calculation, we find that this integral is equal to 



p(u - r', x) In s' + ^ — I - In s' + tp(n + 1) + 7 



71=0 



(u,x" ) 



p=n+l r-\i> ) ) 



(7.25) 



Notice that the first term cancels out the first term to the right-hand side of Eq. (7.24). The rest of Eq. (7.25) is 
simplified by invoking the slow-motion approximation. Because it involves an additional (retarded) time derivative, 
the first term in the sum over p is smaller than p( n+1 > by a factor of order v, and the remaining terms are smaller 
still. Now, Mp(™ +1 ' is already a fac tor of order v 3 smaller than p^ n \ which appears in the direct part of the radiative 
field. This means that in Eq. ( 7.25), the sum over p is 0(v 4 ), and therefore, it will be neglected. 

After substituting Eq. ( 7.25| ) into Eq. (7.24), and expanding the third term to the right-hand side of this equation 
in a Taylor series about n ■ x' — 0, we arrive at 



2M ^ 1 



n=0 



In s' + %p(n + 1) + 7 



(u,x" ) 



+ / p {n+2) («', x' ) In(u - u) du' } (n ■ x' )" dx' . 

J — 00 

This is our final expression for the tail part of the radiative field. 



(7.26) 



D. Total radiative field 



The total radiative field is obtained by adding the direct and tail terms. Combining Eqs. ( 7.18 ) and ( 7.26 ), we 
obtain 



$rad(*,£) 



1 00 1 

-Y- 

r ^— ' nl 

n=0 



/)W(u,f)+2M 



In 



2c 



+ tp(n + l) +7 



p( n+2 \u',x' )du' \{n-x') n dx' 



(7.27) 



We recall that dx' = (1 + 2M/r')d 3 x' , u = t — r — 2Mln(r/c), and that Eq. ( 7.27 ) has been derived on the basis of a 
weak-field, slow-motion approximatio n; it is valid to first order in M, and neglects terms of order w 4 . 

Th e con stant c appearing in Eq. ( 7.27| ) is the same one which enters in the definition of the retarded time it, 
Eq. ( 7.16 ). The radiative field does not actually dep end o n the numerical value of this constant. To see this, let 



(n) 



Ac, where A is a scaling constant. Then Eq. (7.16) implies u 



2M In A, and we have p^ n \u,x') 



{u,x') + 2M In A// ,l+1 ' ) (u, x 1 ) + 0(M 2 ). Substituting these relations into Eq. ( 7.27| ) and discarding all terms of 



P 

order M 2 shows that indeed, <& ra d(cc ) is invariant under this transformation. 

Our current expression for the radiative field has a mathematical structure similar to that of Eqs. (3.11), (3.19), 
and ( 3.20 ), but there appear to be some differences. We now show that these are only apparent, and that in fact, 
Eq. ( 7.27 ) is entirely equivalent to the results of Sec. III. 

We first re-introduce the spherical harmonics, with the relation J41| 



(n-x') n = 4wnlr' n ^ ^ 



2 i(n±l)\ 

v 2 r . , s Y lm (9' <b') Yi m <9, 0), 

( n + l + l )!(a=i)! V ,V ' V ,w ' 

t— rn— — l v ' \ Z J 



(7.28) 



where the sum over I includes even values only if n is even, and odd values only if n is odd. The angles 9' and <fi' are 
the polar angles of the source point a?, and 
reordering the sums, we obtain 



and 4> belong to the field point x. Substituting this into Eq. (7.27) and 
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4tt 



J2 E Yim{e',4>')Y lm {e,4>) £ 



1=0 m=-l 

2M 



P (n) (u,x') 



, u — u 1 

In — — + ip(n + 1) + 7 
2c 



\ (n + I + 1)1(^1)1 

p {n+2 \u',x')du'\r ln dx' 



(7.29) 



The slow-motion approximation now demands that we keep only the first two terms (n = / and n = I + 2) in the sum 
over n. After some algebra, we arrive at 



(7.30) 



which is the same as Eq. (3.11). Here, 

Z lm (u) = Z lm (u) + 2M 

and 

Zi m (u) 



ln[^-^) +r(/-:- 1) - 



Zi m (u')du', 



-in 



(2Z + 1)!! \du 



d 



{rd u f 
2(21 + 3) 



+ 0{v A ) 



p(u, x) r l Y lr , 



i dx. 



(7.31) 



(7.32) 



This is almost, but not quite, the same as Eqs. ( |3~19| ) and ( |3~20| ). 

To properly compare our results with those of Sec. Ill, we must account for the different choices of coordinate 
systems. The coordinates used in this section, and those for which Eq. (7.4) holds, are the harmonic coordinates 
{t, x, y, z}. From these we have constructed the spherical coordinates {t, r, 9, </>} in the usual way, and in this coordinate 
system, dx = (1 + 2M/r)r 2 dr dfl, where dJl = dcos9d<fi. These coordinates are distinct from the Schwarzschild 
coordinates used in Sec. Ill, which we now denote {t, f, 9, <fi}. The transformation between the two coordinate systems 
is M 



t = t, r = r + M, 



(7.33) 



We therefore have dx = r 2 dr df l = dx, which is the volume element of Sec. III. We also have r = r (1 — IM/r), and 



substituting this into Eq. (7.31) yields 



Zi m (u) = 



47T 



(2Z + 1)!! \du 



1+ 2(2iT3)" y + ° iV \ 



p{u,x) r l Yi m (9,ip)dx. 



(7.34) 



This is the same as Eq. ( 3.20|) . 

Finally, a specific choice for c can be made by demanding that u = t — r — 2Mln(r/c) be equal to u = t — f — 
2Mln(f/2M), which is the retarded time encountered in Sec. III. (We have approximated f/2M — 1 by f/2M in the 
logarithm.) A short calculation gives 



c = 2Me- 1/2 , 



and with this choice, Eq. (7.31) becomes 
Zim{u) = Z hn (u) + 2M 



, u — u'\ , . 1 

In + tb(l + 1) H h 7 

AM J V 1 ' 2 ' 



Zi m (u')du'. 



(7.35) 



(7.36) 



This is the same as Eq. ( |3.19 ). We therefore have precise agreement between the results of this section and those of 
Sec. III. 
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APPENDIX A: TENSOR FIELDS ON A TWO-SPHERE 



We gather, for the benefit of the reader, several results pertaining to the "edth" differential operators and the asso- 
ciated spin- weighted spherical harmonics. The discussion follows closely Ref. (43|, but it is essentially self-contained. 
We consider S 2 , a spherical two-dimensional space with metric 

ds 2 = g ab dO a d6 b = r 2 (d0 2 + sin 2 9 dip 2 ), (Al) 

on which fields of various tensorial types are defined. For simplicity, all geometric objects considered in this Appendix 
will be confined to this space. However, extension of the following considerations to four-dimensional, spherically 
symmetric spacetimes is immediate. 

It is useful to introduce a set of basis vectors, m a and rh a (a bar denotes complex conjugation), which satisfy the 
relations 

m a m a — fh a m a = 0, m a m a = 1. (A2) 

A particular choice is 

T T 

m a = — -= (l,isin0), m a = — -= (1, — isin0). (A3) 
v 2 v 2 



This choice is not unique, because the relations (A.2) are preserved under the transformation 

m a -» e** m a , fh a -» e-^fh a , (A4) 

where ip is any constant. 

We may use the basis to decompose tensor fields. For example, a vector field V a may be expressed as 

V a = -iVm a + iVm a> (A5) 

where 

1 V = V a m a , _iV = V a m a . (A6) 
Similarly, a symmetric tensor field T ab is decomposed as 

Tab = -2Tm a m b + 2 Tm (a m b) + 2 Tm a m b , (A7) 

where 

2 T = T ah m a m b , Q T = T ab m a rh bl - 2 T = T ab m a m b . (A8) 



The spin-weight of a field is determined by how the field transforms under (A4). By definition, a field has spin- weight 
s, and is denoted s r], if 

sV - e^.ri (A9) 

under the transformation. For example, -\V has spin-weight s = —1, while 2 T has spin-weight s = 2. 
The covariant derivatives (with respect to g ab ) of the base vectors are given by 

m a , b = — cot 9 m a (m b - fh b ) (A10) 

V2r 

and its complex conjugate. It follows that the covariant derivatives of arbitrary tensor fields can be conveniently 
expressed in terms of the "edth" differential operators, d and d, which are defined by 

/ d d \ ~ ( d d \ 

d = - — +icsc6— -scotO , d=- — -icsc0— + scot0 . (All) 
\o0 dip J \o0 dtp J 

It should be noted that these operators depend on s, the spin-weight of the object on which they act. (The original 
notation Q for these operators was d~ and <9", respectively.) For example, 
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Va: b = (d-iV)m a m b + (d-iV)m a mb + (diV)m a m b + {diV)m a mb ■ (A12) 

From this relation it is clear that d raises the spin- weight by one unit, while d lowers it by one unit. For example, 
d\V = — \/2rV aib m a m h has spin-weight s = 2. 

The "edth" operators can be manipulated efficiently when working under an integral sign. Given two smooth, 
complex functions s —if(0,(f>) and s g(6,(j>), the following identities are easily established by straightforward partial 
integration: 

J@.-if).gdn = -J a -if(Blg)d£l, J (d s g) s .Jdn = - J s g{B^J)<m, (A13) 
where dfl — d cos 6 dip. 

The "edth" operators can be used to generate sets of spin- weighted spherical-harmonic functions, denoted s Yi m (9, 4>). 
Each set (corresponding to a fixed value of s) is complete, and members of a given set obey the usual orthonormality 
relations. The defining relations are oYi m = Yi m (the usual spherical harmonics), and 



8 s Y lm = y/{l-s)(l + s + l) a+1 Y lm , d s Y lm = - v /(/ + s )(/- s + l) s _ 1 y Im . (A14) 

The spin-weighted spherical harmonics also satisfy the relations 

-s%- m = (-l) s+m s Y lm . (A15) 

The spin-weighted spherical harmonics can be combined with basis vectors to form tensorial spherical harmonics 
j28j. For example, 

Y^ lm = -j=(-iY lm m a - x Y lm fh a ), Y a B ' lm = --^(^Y^ m a + x Y lm fh a ), (A16) 

are vectorial spherical harmonics. The superscript E indicates that under a parity transformation, Y^' m has electric- 
type parity: Y^' lra — > (— XyY^' ; the superscript B indicates a magnetic-type parity: Yj 3,1 " 1 — > (—ly^Y^' . 
Similarly, 

T ab" lm = -^=(-2^m m a m b + 2 Yi m fh a fh b ) , T^ 2 ' 1 " 1 = --^=(_ 2 Y5 m m a m b - 2 Yi m fh a fh b ) , (A17) 
are tensorial spherical harmonics. 



APPENDIX B: EVALUATION OF TWO FUNCTIONS 

We evaluate, in the limit z — > 0, the functions Ai(z) and -B;(z) defined by Eqs. (2.14) and ( 2.15| ). 
To evaluate Ai(z) is easy. By using the expansions Si(2z) — 2z + 0(z 3 ) and z l n v j v — —z/(2p + 1) + 0(z 3 ), we 
quickly arrive at 



Ai(z) = z - z 



(Bl) 



The sum evaluates to 1 — 1//, and Ai(z) reduces to the result quoted in the text — Eq. ( 2.16| ). 

To evaluate Bi(z) requires more work. We begin by recalling the series expansions for the cosine integral (Ref. Q 
p. 232), 



Ci(2z) = 7 + ln(2z) + ^ 



(-1)" 
' 2n(2n) 



{2zY n , 



(B2) 



the squared sine, 



sin z = — 



E 



J 2(2n)\ 



(2zf 



(B3) 
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and the squared spherical Bessel functions of the first kind (Ref. pp| , p. 960), 

(-l) fe r(2p + 2 + 2/fc) /z\2 P - 



^fe!r(2p + 2 + A)r»( P +§ + *) 

Substituting these into Eq. ( |2.15 ) gives the series 

oo 



(I)' 



where, after some rearranging, 

{~l) n+1 2 2n {n- 1) 7rr(2n) 



1 1 



2n(2n)! 2 2 ™r 2 (n + §) ^ \P p + 1/ (n -p- l)!r(n + p + 1)' 

with P = min(n — 1,1 — 1). Additional manipulations bring 6„ to the form 



(-l) n+1 2 2n (n- 1) 
2n(2n)! 



(-l)P 



The sum evaluates to 



n - 1 



n — l-^-^yp p + 1 J (n — p — l)!(n + p)\ 

(-l)£(n-P-l) 
n(P + l)(n + P)!(n-P- 1)!' 



which implies that 6„ = for n < I (because P = n— 1), while 

(_l)n+I 2 2n-l(„_ !)|2 



&n = 



/(2n)!(n-Z- l)!(n + Z-l)! 
for «>/ + !. It is then easy to show that Eq. (B5) reduces to the result quoted in the text — Eq. (2.16|). 



(B4) 



(B5) 



(B6) 



(B7) 



(B8) 



(B9) 



APPENDIX C: REGGE- WHEELER FUNCTION IN THE ASYMPTOTIC LIMIT 



We wish to evaluate the Regge- Wheeler function X[*(z), as given by Eqs. ( |2.8| ) a nd (|2.13 ) (with a = and b — —7), 
in the asymptotic limit z — > 00. By com paring with the low- frequency limit of Eq. ( |2.4| ), we will then compute A\ n in 
the normalization prov ided b y Eq . (2 



To express Eqs. ( p.8|) and ( |2.13| ) in the limit z — > 00, we use such asymptotic relations as Si(2z) ~ 7r/2, Ci(2z) ~ 0, 
z 3 (nd p - jin p )j p ~ gpL - (-i) l ~ p ]zn h ~ and ji+i ~ After some algebra, we obtain 



X," ~ (1 - e\)zj x - e[ln(2z) - 0i]zrn + 0(e 2 ), 



where (ii is given by Eq. (2.18). 

We must now compare this result with the low- frequency limit of Eq. (|2.4|), which we rewrite as 



H 



Afe 



- A? ut e lz ' 



(CI) 



(C2) 



where z* = z + eln(z/e — 1). Expanding the phase factors in powers of e, and using the asymptotic relations 



where 



(±i) l+1 (zji ± izni), yields 



Xf - (1 + eAf)zji + s(Ai - \nz)zm + 0(e 2 ), 



1 + eAj = (- 



ieA l = 



l+l j^in gie In e 



l+l /lout — ie\ne 



+ (^) i+ M^ t e 



l+\ j^ui^is In e 



i+1 yiout — ielne 



(^) i +M^ t e 



Finally, comparing Eqs. (ClJ) and flC3|), and using the relations (C4), we arrive at 

AT = l{i) l + 1 e - ie( - in2e -^ [l - f e + 0(e 2 ) . 
From this and Eqs. (§1|, fl2.13| ), we obtain Eq. ( [2.17D . 



(C3) 



(C4) 



(C5) 
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APPENDIX D: INVERSE FOURIER TRANSFORM OF TAIL CORRECTIONS 



We wish to take the inverse Fourier transform of the function 



(Dl) 



where c is a constant and F(u>) an arbitrary, square-integrable function. In other words, we wish to compute the 
function T{u) given by 



f(uS)e 



1 du. 



(D2) 



We shall do so in the spirit of the slow-motion approximation, by formally treating M as a small parameter. We 
follow closely the derivation found in Appendix A o f Ref . [[l9| . 

We first expand the exponential factor in Eq. (Dl) to linear order in M, and combine the result with the (l + 7rM|w|) 
factor. We then substitute the identity iojln \ui\ + tt\u>\/2 = iu>\n(—iu>). After a few lines of algebra, we obtain 



T{u) = F{u) - 2M(ln4M - c)F{u) + 2M / F(lu) iuhx{-iu)e 



where F(u) is the inverse Fourier transform of F(uj), and a dot indicates differentiation with respect to u. 
To evaluate the integral, we write ln(— iuj) in a different form by using the identity (Ref. pQ ], p. 573) 

/>oo 

In \i = —7-/1 / e~^ x h\xdx 1 



(D3) 



(D4) 



with fj, = — iu>; 7 ~ 0.57721 is Euler's number. Strictly speaking, this identity is valid only if the real part of fi is 
positive. This problem can be circumvented by introducing a regulator e > 0, and setting /1 = — iuj + e. The limit 
e — > can be taken after integrating over u>, which yields 

(D5) 



F(u) = F(u) -2M(ln4M - c-j)F(u) + 2M / In x F(u - x) dx. 

Jo 



We write this in its final form as 

T{u) = F(u) + 2M 

This is the desired result. 



In +c + 7 

4M 1 



F{u') du'. 



(D6) 



APPENDIX E: CHANDRASEKHAR TRANSFORMATION FOR THE ELECTROMAGNETIC FIELD 



We deriv e th e relation between i?p(w;r), the solution to the homoge neo us version of Eq. (L2) with bo un dar y 
conditions (4.9), and X/ ? (o;;r), the solution to the s — 1 version of Eq. (2.1) with bounda ry c onditions (2.3) ( |2.4| ). 
For convenience, in this Appendix we set to unity the arbitrary constant appearing in Eq. ( |2.3[ ). 

Direct substitution shows that if Xi(cu;r) satisfies the generalized Regge- Wheeler equation (with s = 1), then 
Ri(uj;r) = r£Xi(cu;r) satisfies the homogeneous Teukolsky equation. Here, C = fd/dr + iu>. The desired relation 
must therefore have the form 



1(1+1) 
1 - AiMui 



1 + o(/ 2 ) 



(El) 



where the constant x must be chosen so that the normalization of Rf 1 (uj; r) agrees with Eq. (4.9). 

To find x an d to relate Q\ a (uj) to A\ n (uj), we need expressions for Xf 1 (ui; r) that are more accurate than Eqs. (2.3) 
and (2.4). By solving the generalized Regge- Wheeler equation, we find that 



(E2) 



near r = 2M , while 
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I f H ( W ;r)=irH 1 + 



1(1 + 1) 
2iujr 



+ 0[(cjr) 



— 21 l_ — iujT* 



(E3) 



near r — oo, where the dots designate terms proportional to e iu ' r . We are now in a position to verify that near 
r = 2M, rCe- l ^ r ' ~ 1(1 + l)fe~ lulr */(l - 4iMu), and that near r = oo, rCe- wr * = + l)(iwr)- x e~ iur " . 

Combining these results with Eqs. (||J|), IPI ; and @)' we nnd 



X 



1 - 4iMw 
/(/ + !) 



(E4) 



and 



(1 -4iMw)^l n (w). 



(E5) 



Equation (4.15) follows immediately. 
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